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Abstract 

We introduce the generalized spatial discretization of the Kardar-Parisi-Zhang (KPZ) equation 
in 1+1 dimensions. We solve exactly the steady state probability density function for the discrete 
heights of the interface, for any discretization scheme. We show that the discretization prescription 
is a consequence of each particular model. From the ballistic deposition model we derive the 
discretization prescription of the corresponding KPZ equation. 
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The interface growth models has attracted much attention during the two last decades 
due to its widespread application to many systems , such as film growth by vapor or 

chemical deposition, bacterial growth, evolution of forest fire fronts, etc. For such systems, 
the major effort has been concentrated in the identification of the scaling regimes and their 
classification into universality classes through Monte Carlo simulation of the discrete models 
or renormalization-group analysis of the continuous equations that describe the evolution of 
the interfaces in the coarse-grained approximation. This powerful tool allows a correspon- 
dence between the discrete models and the continuous equations. Another useful method 
is to derive continuous evolution equations from the transition rules of the discrete growth 
models based on a regularizing scheme and coarse-graining of the discrete Langevin equa- 
tions 0, O, . Phenomenological equations, selected according to symmetry principes and 
laws conservation, are often able to reproduce various experimental data. Between these 
phenomenological equations the introduced by Kardar, Parisi, and Zhang (KPZ) has been 
succesful [7] to describe the properties of rough interfaces, but is also related to Burgers 
equation of turbulence, and to directed polymers in random media [2|. The KPZ equation 
describes the evolution of the profile h(x, t) of the interface at position x and time t: 



Oh d 2 h X (dh\ . . 
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where v and A are the diffusion and nonlinear coefficients, respectively. The Gaussian 
thermal noise rj(x,t) has zero mean and covariance 



(r)(x, t)7](x', t')) = 2e5(x- x')5(t - t') 



(2) 



where e is the noise intensity. Here and elsewhere ( ) denotes average over the noise real- 
izations. The KPZ equation differs from the Edwards- Wilkinson (EW) equation [8], that 
describe interfaces growing under the effect of random deposition and surface tension, in 
a nonlinear term due to microscopic lateral growth. The direct numerical integration is a 
powerful and simple tool to compute the exponents that characterize the universality class 



of a given continuous equation 9|, [ljj, [12j. In this context, the KPZ equation has been 
integrated using a unusual discretization method with exact steady state probability density 
function jl^j], a pseudospectral discretization method and a least-square error method 
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although the theoretical justification is not clear. 



The main goal of this Communication is to introduce the generalized spatial discretization 
of the KPZ equation in 1 + 1 dimensions. We show that exists the steady state probability 
density function of the discrete interface heights compatible with the continuous ones, after 
coarse grained approximation, independent of the discretization scheme. Finally, we find 
the discretization prescription for the KPZ equation derived from the ballistic model. 

Generalized discretization. We propose as generic spatial discretization of the KPZ equa- 
tion the following Langevin equation 



dhi 
~dt 



(7) 
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where hi(t) = h(ia, t) is the interface height at the ith-lattice point (i — 1, . . . , N), a is the 
horizontal lattice spacing, and L = a N is the lattice size. Periodic boundary conditions are 
assumed, i.e. ho = h^. Without loss of generality we take the horizontal and vertical lattice 
spacing to be equal. Introducing in Eq. (jSJ) the addimensional difference of heights 
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where £, k = —1, 0, 1 (£ 7^ k), the standard discretized diffusive term is given by 

Li = \ {H? 1 - HQ . 
Our discretized nonlinear term (with < 7 < 1) in Eq. Q is defined as 



N, 



(7). 
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The noise r\i has zero mean and covariance 
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Expanding Eq. (0) and Eq. (jUJ) around the zth-site, the discretized diffusive term and the 
nonlinear term in the 7-discretization are 
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respectively. Several discretizations will produce different results in the roughness even- 



though the difference in numerical accuracy of the height profile are small 



13|. Most nu- 



merical studies are done with the discrete spatial version of the KPZ equation corresponding 



to the usual choice 7 = 1 in Eq. © called standard or post-point discretization. The nonlin- 
ear term = ^(H^l) 2 only depends on the nearest neighbor sites height and minimizes 
the error in approximating d 2 h/dx 2 [see Eq. (J5])]. Moreover, the choice 7 = in Eq. (ffij). 
called anti-standard or pre-point discretization, used here afterwards, corresponds to the 
arithmetic mean of the squared slopes around any interface site. On the other hand, Lam 
and Shin [3| introduced the spatial discretization corresponding to the choice 7 = 1/2 in 
Eq. © that enables an elegant analytical treatment. However, this choice is unusual and is 
only supported by the existence of a steady state probability density function equal to those 
of the linear case A = 0. We explain this special choice and show, through a general calculus 
of the steady state solution, that the generalized discretization has an unambiguous limit in 
the continuous. 

Steady state density. The main feature of the generalized discretization of the KPZ 
equation is that the corresponding steady state probability density function P(h) of the 
discrete heights h = {hi} exists. This density function give rise to the known steady state 



probability densi 



y functional P[h] of the continuous interface height h(x,t) related to the 



KPZ equation 



a. 



Conversely, P(h) is not a direct consequence of the P[h] as we show 
bellow. The probability density function P(h,t) of the discrete interface evolves according 
to the following Fokker-Planck equation 

N Q 7(7) 
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where Jj (h, t) is the probability current density. Replacing in Eq. the general 

P(h) = exp (-- e £ aU^ , (10) 

where is a "potential" function to be derived bellow, we obtain the steady state prob- 
ability current density function 

^ (7) (h)= (uL i + ^N^ + u^jP. (11) 

Lam and Shin showed that P(h) for the linear case (A = 0) is also solution of the 
nonlinear one in the midpoint discretization. We explain their result setting U- = in 



Eq. (fTTj). i.e. the solution corresponding to A = 0. It is easy to show that, under periodic 
boundary conditions, the current j\ = (A/2) P is conserved only if 7 = 1/2 for A 7^ 
(i.e. Y^iLidJi / dhi = 0). Our goal is to find a steady state solution, independent of 
the choice of the discretization and therefore for any nonlinear coefficient A, with constant 
J^(h). It is easy to show, using Eq. (J3J) and Eq. (jUJ), that 
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Using Eqs. (|T2"j) in Eq. (JTTJl we obtain that 

ut^N^+aL?, (13) 

where 
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Eq. (JT3|) gives a steady state solution of the Fokker-Planck equation [Eq. (fTUjl] with conserved 
current 

A dimensional analysis shows that e/v and z//A are proportional to a, and therefore a is 
proportional to a 3 . Besides, note from Eq. (jSJ) that the errors of iV"/ and Lj in approximating 
(dh/dx) 2 and d 2 h/dx 2 , respectively, are at most proportional to a 2 . Thus, we conclude that 
the error of in approximating (dh/dx) 2 is proportional to a 2 [see Eq. Q13|)]. The midpoint 
solution U^ 2 ^ = N- QS) is symmetric under the interchange -£^ i+1 <-> but the term crL 3 

breaks weakly this symmetry in Eq. (|13|1. Notice that the limits A = (EW equation), 
7 = 1/2 (midpoint discretization), and a = (coarse grained approximation) in Eq. ()13j) 
are equivalent between them. 

If we denote by P([h],t) the probability density functional of the interface position func- 
tion h(x, t) in 1 + 1 dimensions, the corresponding Fokker-Planck equation is 
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from which we can obtain the well-known steady state solution 



P[h] = exp 



v r L (dhS* 
e Jo \dx 



Incidentally, this is the solution in absence of the nonlinearities. Moreover, 

J t"i = i© 2p = CDn8t 

for the nonlinear case. A special feature of our discretization is that P(h) is a steady 
state solution for all values of the A. While the nonlinear coefficient is present in P(h), 
it is absent in P[h], because lim a ^o X)£Li a = J L (dh/dxf dx. This means that the 
coarse grained approximation erase in P[h] our knowledge on the discrete model and the 
nonlinearity. Again, for any discretization, the coarse grained approximation preserves the 
nonlinear dominant term of the steady state current, since lim a ^o — J[h]. 

Ballistic deposition model. The KPZ equation is tractable by direct numerical integration 
when we specify its associate spatial discretization. In order to verify this fact, we derive 
the Langevin equation [Eq. Q] for the ballistic deposition (BD) model. The procedure 
used here is based on regularizing the step functions included in the growth rules of the 
microscopic model, in order to obtain the discrete Langevin equation and, afterward coarse 
graining, the KPZ equation. Let us first introduce the general treatment. The discrete 
interface growth at an average rate r and the interface height at an ith-site increase in 
hi(t+r)-hi(t) = a £f=i rf\ where are the rules of the deposition processes. Expanding 
hi(t + r) up to second order in Taylor series around r, we obtain hi(t + r) —hi(t) ~ rdhi/dt. 
Thus, the height evolution equation of the zth-site is given by the Langevin equation 

^ = *J" + *(*), (M) 

where the Gaussian thermal noise 77, has zero mean and covariance 

{ m {t) m {t')) = K® 5{t-t'). (15) 

The first and second moments of the transition rate, in terms of growth rules, are given by 

KP = - E 4 j) , Kf = a5 ijK ?\ (16) 
T 3=1 

respectively. In the BD model, a particle is released from a randomly chosen lattice position 
i above the interface, located at a distance larger than the maximum height of the interface. 



The incident particle follows a vertical straight trajectory and sticks to the interface at time 
t. The height in the column i is increased by max[/tj_i, hi + 1, hi+x]. For this model the rules 
can be summarized as: 
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where 0(z) is the unit step function defined as 6(z) = 1 for z > and 6(z) =0 for z < 0, 
and 5(^,0) = B(z) + Q(—z) — 1 is the Kronecker delta. The representation of the step 
function can be expanded as Q(z) = Y^k=o c k zk providing that z is smooth. In any discrete 
model there is in principle an infinite number of nonlinearities, but at long wavelengths the 
higher order derivatives can be neglected using scaling arguments, since one expect affine 
interfaces over a long range of scales, and then one is usually concerned with the form of 
the relevant terms. Thus, keeping the expansion of the step function to first order in its 
argument and replacing the expansion in Eq. (fTTj) and (jl6j) the first moment is 
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where 
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The average driving velocity vq can be subtracted in the expression of the first moment given 
by Eq. (fTSj) . choosing adequately a moving reference frame. Retaining only the constant term 
in Eq. (fTHj) we obtain K^' ~ a5ijV = eb^ja. Replacing in Eq. (fTo]) we recover Eq. ((7j) 
with noise intensity e = a 2 1> . Notice that in order to integrate numerically the continuous 
equation, we need a continuous representation of the 6-function to numerically compute the 
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coefficients Co and c\ related to the ones of the KPZ equation through Eq. (JTHJ), such as the 
shifted hyperbolic tangent representation 16]. Our parameter of discretization is related to 
the parameters of the microscopic model and the KPZ coefficients [see Eq. ()19|)]. Choosing 
7, the KPZ coefficients A, v and e depends only of one of the microscopic parameters, e.g. 
with 7 = 1/2 the parameters c and C\ are related as c + C\ = 1/2. 

In summary, we propose a finite difference discretization criterion for the numerical in- 
tegration of the KPZ growth equation starting from the corresponding generalized discrete 
Langevin equation. We show that, for any discretization scheme, exists a steady state proba- 
bility density function of the discrete interface heights compatible with the continuous ones, 
after coarse graining. Besides, we derived the KPZ coefficients of the BD model as an ex- 
ample of our generalized discretization. Finally, our results can be used as a new tool for 
the direct numerical integration of growing continuous equations with nonlinear terms due 
to lateral growth like the KPZ one. 
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